Temperature and time scaling of the peak-effect vortex configuration in FeTe .7Se .3 
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An extensive study of the magnetic properties of FeTeo.7Seo.3 crystals in the superconducting state 
is presented. We show that weak collective pinning, originating from spatial variations of the charge 
carrier mean free path (81 pinning), rules in this superconductor. Our results are compatible with the 
nanoscale phase separation observed on this compound and indicate that in spite of the chemical 
inhomogeneity spatial fluctuations of the critical temperature are not important for pinning. A 
power law dependence of the magnetization vs time, generally interpreted as signature of single 
vortex creep regime, is observed in magnetic fields up to 8 T. For magnetic fields applied along 
the c axis of the crystal the magnetization curves exhibit a clear peak effect whose position shifts 
when varying the temperature, following the same dependence as observed in YBa2Cus07_i. The 
time and temperature dependence of the peak position has been investigated. We observe that the 
occurrence of the peak at a given magnetic field determines a specific vortex configuration that is 
independent on the temperature. This result indicates that the influence of the temperature on the 
vortex-vortex and vortex-defect interactions leading to the peak effect in FeTeo.7Seo.3 is negligible 
in the explored range of temperatures. 

PACS numbers: 74.70.Xa, 74.25.Ha, 74.25.Wx 



I. INTRODUCTION 

The study of the vortex properties in type-II supercon- 
ductors is of extreme interest both for investigating the 
basic physics of the superconductivity and for evaluating 
the quality of the materials in view of practical appli- 
cations. Many aspects of this subject are still topical, 
especially with regard to the high temperature supercon- 
ductors (HTS) and the recently discovered iron-based su- 
perconductors. To date, five families of Fe-based super- 
conductors have been discovered: REOFeAs, ("1111", 
RE=r&re earth),i ^Fe 2 As 2 ("122", ^=alkaline earth), 2 
XFeAs ("111", A=Li; Na)£i Fe(Se,Ch) ("11", Ch=S, 
Te)^ and the most recently discovered "21311" fam- 
ily of Sr 2 M0 3 FePn (M=Sc, V, Cr and Pn=pnictogen),I 
Among these families, iron chalcogenides are considered 
of particular interest because of their simple crystal struc- 
ture consisting of Fe ions tetrahedrally coordinated by 
Se and Te arranged in layers stacked along the c-axis, 
without any other interlayer cations, as occurs in the 
pnictides. For this reason, iron chalcogenides are gen- 
erally considered an ideal candidate for understanding 
some open issues of high-temperature superconductivity. 

One of the most intriguing phenomena observed in the 
study of the vortex properties in type-II superconduc- 
tors is the so called peak (or fishtail) effect. A second 
peak in the magnetization curve M(H) of low-T c super- 
conductors has been observed as early as in the sixties 8 
and a first theoretical explanation has been suggested by 
Pippard in 1969.— The peak has been found in proxim- 
ity of H c2 in Nb 3 Sn^ CeRu 2 ^ NbSe 2f ii VsSii 2 . and 
MgB 2 3iii 



The peak effect has been clearly observed also in the 
cuprate superconductors but its features vary from 
one compound to the other. In YBa 2 Cu 3 07_5 crys- 
tals, the peak is broad and temperature dependent; 
its origin has been attributed to the presence of inho- 
mogeneities in the oxygen content^— and associated 
to different mechanisms J^r— In the more anisotropic 
Bi 2 Sr 2 Ca„_iCu rl 2rl +4+d the peak in the magnetization 
loop has usually been associated to a 3D to 2D tran- 
sition of the vortex structure! 23 ' 24 However, significant 
differences in the second-peak features have been ob- 
served within this family: in the highly anisotropic Bi- 
2212 and Bi-2223 the peak is sharp and temperature 
independent 2 ^— while in Bi-2201 crystals it exhibits a 
strong temperature dependence^ Regarding the iron- 
based superconductors, the peak effect has been observed 
in FeAs- 1111 ^rM, in FeAs- WM^. and in FeTei _ x Se^ ^ 
even if a clear understanding of the physical origin of the 
phenomenon in this family has not been achieved, yet. 
However, since it has been shown by small angle neutron 
scattering j 4 ^— Bitter decoration 4 ^ and magnetic force 
microscopy 4 ^ that an ordered vortex lattice is not present 
in these compounds, the origin of the peak effect due to 
a phase transition of the vortex lattice seems to be ex- 
cluded. 

Beyond the investigation of the magnetization curves, 
an useful approach for investigating the vortex dynamics 
is the study of the relaxation processes of the magne- 
tization. The critical state of the vortex lattice, which 
determines the hysteresis of the magnetization in type-II 
superconductors, is a metastable state. It follows that 
vortices tend to hop out of their pinning-potential well 
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in order to reach the configuration of absolute minimum 
energy. Such motion usually arises from thermal activa- 
tion, but it can also arise from quantum tunneling (at low 
temperatures) or can be stimulated by external perturba- 
tions, such as microwave shaking of the vortex latticed 
Magnetic relaxation processes have been observed in var- 
ious low-temperature superconductors^ 4 - However, the 
subject has become of even greater interest after the dis- 
covery of high-T c superconductors, because of the higher 
operating temperatures and of the small activation en- 
ergies related to the short coherence length and large 
anisotropy. The concept of thermally-induced hopping 
of the flux lines has been first treated by Anderson and 
Kim i 45 ' 46 In the framework of their model a logarithmic 
dependence of the magnetization (M) on the time (t) 
is expected. This behavior has been verified in various 
superconductors, both low T c and high T c ; however, in 
many experiments, deviations from the logarithmic de- 
pendence of M have been observed, indicating a failure 
of the approximations that bring this expectation^ 4 -^ 
In the framework of the collective pinning theory a more 
complex expression for the dependence of the potential 
energy barrier height (the so-called "interpolation for- 
mula") is proposed: U{J) = Uq/ f-t[(J c / JY — 1] where /x 
is a parameter varying as a function of the vortex- vortex 
and vortex-pinning center interaction! 4 ^ 49 For example, 
in three dimensional systems /i = 1/7 when the creep is 
dominated by the motion of individual flux lines, /i = 3/2 
in the case of collective creep of small bundles and [i — 
7/9 when the bundle size is much larger than the Larkin 
correlated volume^ 9 -. In 1991, Vinokur, Feigel'man and 
Geshkenbein proposed a theoretical model for the ther- 
mally activated flux creep, assuming a logarithmic de- 
pendence of the activation energy on the critical current, 
U(J) = U a ln(J c / 'J). 50 This dependence of U over J is 
a good approximation for the creep activation barrier in 
the single vortex creep regime (limit [i — > 0) and provides 
a proper fit of the experimental U( J) dependence ob- 
served in different superconductors The divergence 
of U (J) when J — > can be understood in the context of 
the collective pinning considering that in the vortex-glass 
state vortex motion is possible only in the presence of a 
current. 44 In the case of this logarithmic dependence of 
the activation energy on the critical current, a power-law 
time dependence of the magnetization is expected^, 

In this paper we study the magnetic properties of 
FeTeo.7Seo.3 crystals in the superconducting state by in- 
vestigating the dependence of the magnetic moment (m) 
on the applied magnetic field (H) the temperature (T) 
and the relaxation time (t). In Sec. [IT] we report details 
on the examined samples and on the experimental tech- 
niques used. The experimental results are shown and dis- 
cussed in two different sections. In Sec. IIII A[ dedicated 
to the vortex pinning properties, we analyze the m(H) 
curves and deduce the pinning mechanism in the inves- 
tigated crystals. The magnetic relaxation processes are 
studied in Sec. IIII Bl Conclusions are reported in Sec. lIVI 
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FIG. 1. Temperature dependence of the DC volume suscep- 
tibility (demagnetization corrected) for the investigated sam- 
ples. 



II. SAMPLE AND EXPERIMENTAL 

The FeTeo.7Seo.3 crystals investigated were grown 
by a modified Bridgman-Stockbarger method, starting 
from a nominal Fe:(Te,Se) ratio of 0.9:1. Details of 
the procedure used for preparing the samples are re- 
ported in Ref. 5 -. The refined composition of the sam- 
ples is Fe1.013Teo.68Seo.32, indicating a very low iron 
excess^ It has been shown that reducing the Fe excess in 
Fei+^Tei-ySey favors the occurrence of superconductiv- 
ity and weakens the antiferromagnetic order.— >5£- How- 
ever, according to the phase diagram of the compound, 57 
a little excess of Fe is needed for stabilizing the structure. 
It is also worth noting that, due to Se-doping for Te, less 
Fe is allowed to occupy the additional site, since both the 
effects of reducing x and increasing y result in shrinking 
and re-shaping the FeTe4 tetrahedra; so, the lower is the 
Se content, the more difficult is to limit the Fe excess^ 

Two different crystals, weighting « 50.5 mg (Sample 
A) and s» 0.7 mg (Sample B), have been investigated in 
the present study. The susceptibility curve x(T) of the 
two samples has been deduced from a zero field cooled 
m{T) measurement performed at 10 Oe by a Supercon- 
ducting Quantum Interference Device (SQUID) with the 
field applied along the c axis. Both curves have been cor- 
rected for the demagnetization effects according to the 
formula reported in Ref.— . The results are shown in 
FigUJ the x(T) curve of Sample A exhibits a supercon- 
ducting transition with T c = T(90%) « 10.6 K, AT C = 
T(90%) - T(10%) « 0.8 K whilst for Sample B it results 
T c = T(90%) w 10.9 K, AT C = T(90%) - T(10%) 
1.0 K. Both the samples exhibit bulk superconductivity 
with x(4.25iT) « -1. 

The samples have been characterized also by x-ray 
diffraction (XRD) both in a powder diffractometer (on 
manually ground crystals) and in a 4-circle diffractome- 
ter. The rocking curves indicate the presence of a dis- 
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tribution of the c-axis crystallites with a full width at 
half maximum (FWHM) of about 1.4° for the Sample 
A and about 0.3° for Sample B. Our results indicate 
that, whereas the sample B is of high crystalline qual- 
ity, the sample A is more likely to be composed of few 
single-crystalline domains well aligned along the crystal- 
lographic c-axis. 

The magnetization curves m(H) of the 
FeTeo.7Seo.3 samples were acquired by means of a 
SQUID magnetometer and a VSM (Vibrating Sample 
Magnetometer). The maximum applied magnetic field 
was 4.5 T for the SQUID and 8.8 T for the VSM. Mea- 
surements were performed in the range of temperatures 
from 4.25 K to T c , with the magnetic field applied both 
parallel and perpendicularly to the c axis. The VSM 
measurements were executed with different values of 
the magnetic field sweep rate (SR) namely 0.05 T/min, 
1 T/min and 2 T/min. Furthermore, by means of 
the VSM, we studied the relaxation of magnetization 
over periods of time up to 6000 sec, for fixed values of 
temperature and magnetic field. In order to perform 
measurements in the trapped flux configuration, the 
relaxation measurements were performed with the 
following procedure: i) the sample is zero-field cooled; ii) 
the magnetic field is increased from zero to 8.8 T with 
SR= 2T/min; iii) the field is decreased from 8.8 T to 
8 T with SR= 0.1 T/min; iv) the magnetic relaxation at 
8 T is recorded for 6000 sec; iv) the field is decreased by 
1 T with SR= 0.1 T/min; v) the magnetic relaxation is 
recorded for 6000 sec. Steps iv) and v) are repeated for 
registering the relaxation down to T, in steps of 1 T. 



III. RESULTS AND DISCUSSION 
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FIG. 2. Magnetic moment vs. external field curves obtained 
at different temperatures in Sample A for H \\ c (a) and H _L c 
(b). The insets show the results obtained in Sample B at 
T = 4.25 K and T = 7 K. 



A. Vortex pinning properties 

Fig. [2] shows the magnetic- moment curves of the Sam- 
ple A acquired with the VSM at different temperatures, 
for H || c (a) and H _L c (b), for a sweep rate of 2 T/min. 
Increasing the temperature, the pinning becomes weaker 
and consequently the width of the hysteresis loop de- 
creases. A second peak in the M(H) curve is clearly 
observed for H || c, up to temperatures close to T c . A 
second peak is also present for H _L c but less pronounced 
and detectable only for temperatures higher than 5 K. In 
the inset of Fig. [2] the results obtained for the Sample B 
at T = 4.2 K and T = 7 K, are reported for the two dif- 
ferent orientations of the field. In this second peak 
in the m(H) curve is not observed for H _L c, whereas 
it is clearly evident for H \\ c, in analogy to what ob- 
served in YBa2Cu307_5 crystals^ and in other Fe-based 
superconductors i 35 ' 59 Since Sample A presents a wider c- 
axis distribution of single-crystalline domains than Sam- 
ple B, we argue that the weak peak effect revealed for 
H _L c in this sample is due to the slight misalignment of 
the crystallites forming the sample. We would also like 
to specify that the experimental m(H) curves of Sample 



B for H || c are qualitatively similar with those of Sam- 
ple A at comparable temperatures. This suggests that a 
misalignment of the layers forming the crystal is crucial 
only for magnetic fields applied perpendicularly to the c 
axis of the sample. The m{H) curves obtained for H _L c 
exhibit a paramagnetic-like background arising from the 
experimental setup. This background is not observed in 
the m(H) curves measured for H _L c in the SQUID. 
In any case, the background does not alter the results 
obtained in the following, which rely only on the sepa- 
ration between the positive and negative branches of the 
magnetic-moment loop. 

From the m(H) curves measured for H \\ c at dif- 
ferent T, we determined the temperature dependence 
of the magnetic field value corresponding to the second 
peak in magnetization (H pea k). The results for samples 
A and B are shown in Fig. [3] It follows that for the 
FeTeo.7Seo.3 superconductor the position of the second- 
peak shifts towards lower fields monotonically on increas- 
ing the temperature, in analogy to what observed in the 
YBa 2 Cu307_,5 superconductor^ For both samples in- 
vestigated, the peak becomes undetectable at T w 10 K. 
The continuous lines in Fig. [3] are the best fit curves ob- 
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FIG. 3. Temperature dependence of the magnetic field value 
at which the peak effect occurs, for the two measured samples. 
Continuous lines are the best fit curves obtained supposing 
the same T dependence observed in YBa2Cu307_a. 
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FIG. 4. Field dependence of the critical current density at 
different temperatures and for H \\ c, in Sample A. The inset 
shows the results obtained at T = 4.25 K and for H \\ c for 
the Sample B. 



tained supposing the same T-dependence of H pea k as ob- 
served in YBa 2 Cu 3 7 _ 5 ^ H peak = A • (1 - T/T*) 3 / 2 , 
where A is a constant and T* is the temperature at 
which the peak is undetectable. The best fit parameters 
(with the relative statistical errors) are: A — 5.9 ±0.1 T, 
T* = 10.3 ± 0.1 K for Sample A and A = 6.0 ± 0.2 T, 
T* = 10.8 ± 0.2 K for Sample B. This remarkable simi- 
larity between the results obtained in the YBa2Cu30 7 _,5 
and in the FeTeo. 7 Seo.3 superconductors, might suggest 
that the peak effect in the two systems has an analogous 
origin. 

The magnetic field dependence of the critical current 
density (J c ) has been extracted from the m(H) curves, 
for different values of the temperature, using the Bean 
critical state formulae . 61 ' 62 For slab samples in a perpen- 
dicular magnetic field: J C (T, H) = 3Am(T, H)/w 2 d(3l - 
w), where Am(T, H) is the separation between the two 
branches of the magnetic-moment loop, I and w are the 
length and the width of the sample (I > w) and d is the 
thickness. The J C (H) curves obtained for the Sample A 
at different temperatures are presented in Fig. 01 in the 
inset of the same figure the J C (H) curve at T — 4.25 K of 
the Sample B is shown. For both samples investigated the 
J c values at 4.25 K are of the order of 10 4 A/cm 2 . This 
value is smaller than what measured in Fei+^Tei-ySe^ 
crystals with higher y values than our samplesiS&S 4 How- 
ever, the J c values become comparable if related to the 
same reduced temperature T/T c £^ The fact that the J c 
values obtained for Sample A are very similar to those ob- 
tained for the high-quality single crystal Sample B, where 
no weak links are present, indicates that weak links ef- 
fects are not important for Sample A either. 

The critical current density J c is always limited by 
the depairing current density Jo = 4£? c /3\/6£toA, where 
B c and A are the thermodynamic critical field and the 
GL penetration depth, respectively^ Important infor- 



mation can be deduced from the ratio J c /Jq. For low- 
T c superconductors pinning is usually strong (J c /Jo ~ 
10~ 2 — 10 _1 ) resulting from the interaction of vortices 
with extended defects such as, e.g., precipitates or grain 
boundaries* ~ On the contrary, for the cuprate family pin- 
ning is usually weak (J c /Jq ~ 10~ 3 — 10~ 2 ) normally aris- 
ing from point defects, e.g., oxygen vacancies.— For the 
FeTeo.7Seo.3 we estimate J c /Jq ~ 10~ 3 indicating that 
also in the case of the Fe-based 11 family pinning is weak. 
The microscopic origin of pinning in Fei+^Tei-^Se^ will 
be discussed in the following. 

As already mentioned, the presence of the peak effect 
in the magnetization curve has been widely documented 
in the literature, both in low and high T c superconduc- 
tors. Its origin has been associated to different processes, 
depending on the particular system investigated^ ' 18 ' 22 ' 29 
In order to shed light on the mechanisms that rule pin- 
ning in the Fei+^Tei-ySej, superconductor, we have in- 
vestigated the magnetic-field dependence of the pinning 
force density, F p . It has been shown for a large variety 
of low-T c and high-T c superconductors that the curves 
of F p vs. H obtained at different temperatures may be 
scaled into a unique curve if they are plotted as a func- 
tion of the reduced field, h = H / Hi rr ^— The empirical 
formula that accounts for the scaling is: 

F p = C * h p (l - h) q , (1) 

where C is a proportionality constant, p and q are two 
parameters whose values depend on the origin of the pin- 
ning mechanism^ The different contributions to flux- 
pinning are usually catalogued into two main categories: 
i) 51 (or normal) pinning, arising from spatial varia- 
tions in the charge carrier mean free path (I); ii) 5T C 
(or 8k) pinning, associated with spatial variations of the 
Ginzburg parameter (fc) due to fluctuations in the tran- 
sition temperature T e i 49 ' 67 In the last category fall also 



non-superconducting metallic particles whose dimension 
is smaller than the coherence length (£) since, in this 
case, the superconductivity is induced by the proximity 
effect. A classification is also made for pinning centers, 
as a function of the number of dimensions that are large 
with respect to the inter- vortex distance d ss (4>o/B) - 5 . 
Following the definition given by Dew-Hughes in Refi^ 7 -, 
in this paragraph we refer to point pins as regions whose 
dimensions in all directions are less than d, line pins have 
one dimension larger than d, grain- and twin- boundaries, 
which have two dimensions greater than d, act as surface 
pins, volume pins have all dimensions large with respect 
to d£L 

In the framework of the Dew-Hughes model^ 7 - differ- 
ent values for p and q in Eq.[I]are expected, as a function 
of the specific pinning mechanism involved. Correspond- 
ingly, the theoretical F p vs. h curves present a maximum 
at different h values. In the case of 51 pinning, the max- 
imum is expected at h = 0.33 [p = 1, q = 2) for point 
pins and at ft, = 0.2 (p = 1/2, q = 2) for surface pins, 
such as grain boundaries (the same dependence has been 
predicted for shear-breaking in the case of a set of planar 
pins££); no maximum is expected in the case of 51 volume 
pinning, being p = and q — 2. The maximum of the 
F p (h) curve is expected at higher h values in the case of 
5T C pinning; in particular it occurs at h = 0.67 (p = 2, 
q = 1) for point pins, at h — 0.6 (p — 3/2, q = 1) for sur- 
face pins and at h = 0.5 (p = 1, q = 1) for volume pins. 
Therefore, important information on the physical origin 
of the pinning mechanisms can be achieved by analyzing 
the scaled F p (h) curves. 

The comparison of the predicted pinning functions 
with the experiments requires an estimation of Hi rr , de- 
fined as the -ff-value at which J c — 0. Starting from the 
experimental J C (H) curves, Hi rr is usually determined 
as the extrapolated zero value in the so-called Kramer 
plot^a where (J c 1/2 • H 1 ^) is plotted as a function of 
H. This procedure has been successfully used in the case 
of wires, poly crystalline samples and crystals j 7 -^— How- 
ever, if pinning is ruled only by defects whose dimensions 
are smaller than the intervortex distance, it is expected 
Jc^ 2 oc (1 — h). 67 In this case, Hi„ may be determined 

1/2 

by doing a linear extrapolation down to zero of the J c 
vs. H curve i 67 ' 74 ' 75 In Fig. [5] (a) and (b) we report the 

1/2 

Kramer plots along with the J c vs. H curves, obtained 
at two different temperatures, for Sample A and Sam- 
ple B, respectively. We observe that the Kramer plot 
presents a wide linear behavior for both samples. On the 

1/2 

other hand, the J c curve exhibits a linear behavior in a 
wide range of magnetic fields, only for Sample B. For this 
reason, we have extracted Hi rr from the Kramer plot for 
Sample A, whereas for Sample B we cannot discriminate 
a-priori which is the most correct procedure. However, 
a more conclusive result can be achieved if one considers 
simultaneously both the J c vs. H dependence and the 
scaled F p (h) curves. If one hypothesizes that pinning is 
ruled by point pins and consequently extracts Hi rr from 
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FIG. 5. Kramer plot (right axis) and J^ 2 vs. H (left axis) 
curves obtained at T = 6 K and T — 8 K for sample A (a) 
and sample B (b), as described in the text. 

1/2 

the J c vs. H curve then the corresponding F p (h) curve 
should exhibit a maximum at h « 0.33. Fig. [S] (a) and 
(b) present the normalized pinning-force density for sam- 
ples A and B at different temperatures as a function of 
the reduced field based on Hi rr values deduced from the 
Kramer plots, whereas the inset of Fig. 6 (b) is based 

1/2 

on Hi rr values determined from a linear J c vs. H de- 
pendence. For the sake of clearness, the superscript of 
and Hf rr in the figure indicate respectively the H irr 
values extracted from the Kramer plot or hypothesizing 
J c 1/2 oc (l-h). 

The pinning force curves of Sample A, obtained at dif- 
ferent temperatures, scale well and present a maximum 
at h w 0.27 while, for Sample B, the maximum occurs at 
h ss 0.3 if one considers H irr = or at h ~ 0.33 if 

H irr = Hfrr- These results indicate that, for the single 
crystal Sample B, the pinning mechanisms is ruled by de- 
fects whose dimensions are smaller that the intervortex 
distance in the investigated field range and as a conse- 
quence, the most appropriate procedure for determining 

1 /2 

Hirr{T) is from the J c vs. H curve. On the con- 
trary, for Sample A the position of the maximum in the 
Fp(h) curve suggests that a contribution to pinning com- 
ing from surface pinning cannot be excluded. From the 
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pinning in the single vortex regime it is expected that 
the critical-current density variation with respect to the 
reduced temperature (r = T/T c ) is described by the fol- 
lowing expression;^ 



FIG. 6. Normalized pinning force density curves as a function 
of the reduced field, H/H^, for sample A (a) and B (b). 

indicates that the irreversibilty field has been deduced 
by a linear extrapolation in the Kramer Plot. The inset in 
Fig. [6] (b) shows the results obtained in sample B, scaled by 
using the irreversibility field, Hf rr , deduced by the jV 2 vs. 
H curves 



scaled pinning-force density curves we can also deduce 
information on the type of pinning centers. In particu- 
lar, the fact that the maximum in the F p (h) curve occurs 
for values of the reduced field < 0.33 is an indication that 
the pinning centers in FeTeo.7Seo.3 are of 51 type. It is 
worth noting that this conclusion does not depend on the 
specific procedure used for determining H/ rr . In fact, in 
the case of pinning contribution coming from 5T c -type 
pins, the maximum is expected to occur at h 0.5^ 
which is far away from what observed in our samples. 

The analysis of the scaled pinning force curves in the 
frame of the Dew-Hughes model revealed to give impor- 
tant information about pinning in low-T c , cuprates as 
well as Fe-based superconductors i 17 ' 34 i 68 ' 76 i 77 . However 
since in our samples pinning is weak ( J c / Jq ~ 10~ 3 ) 
and thermal fluctuations are large (Ginzburg number 
~ 10 ~ 3 ), more conclusive results about pinning origin 
in FeTeo.7Seo.3 are expected to be achieved by analyzing 
the experimental results in the framework of the collec- 
tive pinning theory^ Following the theoretical approach 
proposed by Griessen et al^ in the case of 5/-type weak 



J c (r)/J c (0) = (l-r 2 ) 5 / 2 (l + r 2 )- 1 /^ 



(2) 



while, for 5T C pinning, it is: 

J c (r)/J c (0) = (1 - r 2 ) 7 / 6 (l + r 2 ) 5 / 6 . (3) 

In Fig. we plot the normalized J c (t) data obtained 
at hqH = 0.5 T and /iqH = 1 T for H || c, along with 
the theoretical curves expected within the scenario of 51 
and 5T C pinning. The J c (t) values have been extracted 
from the J C (H) curves obtained at various temperatures. 
Analogous results have been obtained for sample B. A re- 
markably good agreement between the experimental re- 
sults and the 51 pinning theoretical curve is obtained, 
confirming that pinning in the FeTeo.7Seo.3 samples orig- 
inates form spatial variation of the mean free path. 

The weak collective pinning theory considers that pin- 
ning originates from fluctuation in the density and force 
of defects whose dimension is smaller than the coherence 
length. In this case small values for the critical current 
density are expected (J c /Jo ~ 10~ 3 — 10~ 2 ).— On the 
contrary, if pinning was mainly ruled by extended de- 
fects such as, e.g., precipitates, twinning or grain bound- 
aries, or columnar defects such as dislocation lines, higher 
values for the J c / Jo ratio would be expected (J c /Jo ~ 
10~ 2 — 10~ 1 )M A phase separation in Fei+^Tei-ySej, 
single crystals with different x and y values has been ex- 
perimentally observed by scanning tunneling microscopy 
and electron energy-loss spectroscopy 78 i 79 In particu- 
lar, fluctuations of the local Te concentration (even by 
20% from the average local composition) in nanomet- 
ric regions have been observed in FeTeo.7Seo.3 crystals^ 
These fluctuations could be the origin of the observed 
weak pinning. 

In spite of the inhomogeneous chemical distribution, 
our results indicate that local fluctuations of the crit- 
ical temperature are not important for pinning. This 
interesting result is in agreement with tunneling spec- 
troscopy experiments performed at T — 80 K in optimal 
doped FeTeo.55Seo.45 crystals, indicating that the com- 
pound is chemically inhomogeneous but electronically 
homogeneous i 7 - 9 - 

It is worth mentioning that a coexistence of supercon- 
ducting and magnetic orders has been proposed for the 
Fei+^Tei-ySej, system. In particular such coexistence 
has been inferred by Khasanov et al. based on muon-spin 
rotation experiments. 80 On the other hand, Li et al. have 
ascribed the apparent coexistence of superconductivity 
and magnetism in the Fei+2,Tei_j,Se y system to a phase 
separation in the real spaced As already mentioned, a 
phase separation has also been observed in FeTeo.55Seo.45 
crystals by STMi 7 - 9 - Other authors have assumed a spin- 
glass state to exist between the long range antiferromag- 
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FIG. 7. Normalized critical current density data, as a function 
of the reduced temperature (t), obtained at fioH = 0.5 T and 
HoH = 1 T for H \\ c. The continuous lines are the theoretical 
curves expected in the case of SI and ST C pinning. 



FIG. 8. Plot of the magnetization vs. time (on a log-log 
scale) at T — 4.25 K for different values of the magnetic field. 
The inset shows the field dependence of the relaxation rate 
S= |d logM/d logt[. 



netic order (at low Se content) and the bulk supercon- 
ducting state, at high Se content^ The low-field DC 
susceptibility measurements of our samples, reported in 
FigJTJ do not show any anomaly or feature attributable 
to the coexistence of superconductivity and magnetic or- 
der. In a previous article we have deeply investigated the 
effect of the actual chemical composition and its effect on 
the crystal chemistry and on the magnetic and supercon- 
ducting phase diagram.— We have proved the combined 
effect of Fe-excess (x) and Se-substitutions (y) on the 
transition from a superconducting to an antiferromag- 
netic state and highlighted the importance of controlling 
the real composition in a three dimensional phase dia- 
gram. By keeping under control a low Fe-excess, we have 
obtained bulk superconductivity at Se-content even lower 
than 0.3i££ In 2010, Bendele et al£2- has confirmed our 
study on the excess Fe and merged data from previous 
publications in a three dimensional phase diagram, still 
claiming a coexistence of superconductivity and mag- 
netism. On the basis of our experimental results and the 
contradicting results reported in the literature, it is not 
possible to discriminate whether such coexistence occurs 
intrinsically at the atomic level in Fei+zTei-ySe^ or a 
phase separation and local composition fluctuations are 
responsible for the magnetic behavior observed by some 
authors. Further studies on very clean and chemically 
homogeneous samples, carried out with different exper- 
imental techniques, are needed in order to clarify this 
point. 

At the end of our analysis on the pinning properties we 
would like to remind that both the Dew-Hughes and the 
Griessen models have been developed in the single- vortex 
pinning regime, i.e. neglecting inter-vortex interactions. 
The validity of this approximation in our case is con- 
firmed by the magnetic relaxation study reported in the 
next section. 



B. Vortex dynamics properties 

In this section we report a study of the dynamical 
properties of vortices in the FeTeo.7Seo.3 superconductor 
performed by magnetic relaxation measurements. The 
analysis has been limited to Sample A, since the mea- 
sured magnetic moment of Sample B was too small for 
achieving a good resolution in a wide range of times and 
temperatures. Fig. [8] shows the relaxation of the magne- 
tization (M = m/V) normalized to its maximum value, 
obtained at T = 4.25 K for H || c, plotted in a log-log 
graph. At all the fields investigated and for times greater 
than about 50 sec, a linear dependence of M has been 
observed, indicating a power-law dependence of M vs. t. 
An analogous behavior has also been detected at higher 
temperatures and for H _L c. It has been shown that 
deviations from the expected M(t) curves at short times 
could be due to the magnetic field overshoot occurring 
when the external field ramp is stopped. 47 This overshoot 
produces a shielded flux zone in proximity of the surface 
of the sample, which affects the initial relaxation process. 
The inset of Fig. [5] shows the S = |d logM/d log£| values 
obtained at different H values; large values of S have al- 
ready been observed in the Fei+ x Tei_ySe H system as well 
as in other Fe-based superconductors . 35 ' i 85 Since the S 
value is related to the pinning potential energy barrier 
height,— the minimum observed in the S(H) curve is as- 
sociated to the presence of the peak effect in the J C (B) 
curves. 

As already described in the introduction, in the frame- 
work of the Anderson and Kim theory, a linear depen- 
dence of M on log(i) is expected; this result stems from 
two basic assumption: i) the pinning potential energy 
barrier height decreases linearly with the current den- 
sity: U = U Q (1 - J/J c ); ii) U /k b T > 1, which allows 
to hypothesize that the thermal-induced hopping rate is 
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proportional to the Arrhenius factor e^^ u °/ kbT \ Sup- 
posing a linear dependence of U on J is only a first- 
order approximation whose validity has been demon- 
strated to fail many times4^ As reported by Vinokur, 
Feigel'man and Geshkenbein, a power-law dependence 
of M vs. t is expected if a logarithmic dependence of 
the activation energy on the current density is supposed: 
U(J) — Uq ln( J c j J)££ In this case, the following expres- 
sions for the time dependence of M are predicted: 

ln(M) =cost- [xf(t)/d](k b T/U )ln(t/To) ,ioi t < t* 

(4) 

ln(M) = cost - (k b T/U a ) 1h(*/to) , for i > ** (5) 

where d is the thickness of the sample, Xf(t) is the po- 
sition of the flux front and t* is the time at which the 
sample is fully penetratedj 50 i 54 

In order to verify that a logarithmic dependence of 
U(J) is actually present for the FeTeo.7Seo.3 crystal in- 
vestigated, we used the method proposed by Maley et 
al.r^ which allows the determination of the U(J) curve 
from the experimental data obtained at different tem- 
peratures. In their paper, the authors show that choos- 
ing a proper value for a time-independent constant (A) 
it is possible to deduce the J dependence of U (or 
equivalently the U(M) dependence) by plotting U — 
-k b T[\n \dM/dt\ - A] vs. M - M eq , where M eq is the 
magnetization at the equilibrium. This procedure is valid 
under the assumption that the temperature dependence 
of U is weak and that the principal effect of increasing 
the temperature is to produce monotonically decreasing 
initial values of M, which is usually valid for T < T c /2^ 

Fig. [9] shows the curves obtained at fi H — 3 T, 
scaled considering A=28, along with the best fit curve 
deduced supposing the following dependence of U on M: 
U = Uo\n(a/(M — &)); the best fit parameters (with 
the relative statistical errors) are: Uq = 113 ± 4 K, 
a = 5.0 ± 0.2 emu/cm 3 and b = -0.39 ± 0.03 emu/cm 3 . 
Similar values for Uq have been obtained in other Fe- 
based superconductors, indicating that the weakness of 
pinning is a general characteristic of these compounds. 37 
The negative value of b is due to the diamagnetic con- 
tribution of the sample holder. A very good scaling is 
obtained for T < 5.5 K, while at T = 6 K a discrepancy 
between the experimental data and the theoretical curve 
starts to be present. As already mentioned in the intro- 
duction, the logarithmic dependence of U on J is a good 
approximation for the creep activation barrier in the sin- 
gle vortex creep regime. As a consequence, our results in- 
dicate that in the FeTeo.ySeo.3 superconductor the vortex 
motion develops in the single- vortex pinning limit even in 
magnetic fields up to 8 T. This is in general unexpected. 
At high fields, the intervortex distance becomes small 
compared to the magnetic field penetration depth and 
thus one would expect that vortex-vortex interactions 
become important and pinning involves vortex bundles. 
However, single vortex pinning up to 9 T has already 
been observed by Inosov et al. in BaFe2- x Co^As2 sin- 
gle crystals by magnetization measurements, small-angle 
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FIG. 9. Magnetization dependence of the pinning potential 
energy barrier height calculated in the frame of the Maley 
model, scaling the data at different temperatures as described 
in the text. 



neutron scattering, and magnetic force microscopy^ Our 
results suggests that also in the Fei +x Tei„ y Se y system 
the vortex-defect interaction dominates up to high fields. 
This experimental evidence further confirms the valid- 
ity of our study on the pinning properties as reported 
in Sec. IIII Al carried out in the framework of the Dew- 
Hughes model where flux-lattice elasticity effects are ne- 
glected. 

In Sec IIII Al we have shown that the m{H) curves ob- 
tained for H || c are characterized by a second peak whose 
position changes with the temperature. In order to inves- 
tigate the magnetic relaxation effects on the peak effect, 
we measured the field dependence of the magnetic mo- 
ment, by means of the VSM, for different sweep rates 
(SR) of the magnetic field. It has been demonstrated 
that the dependence of the hysteresis amplitude on the 
sweep rate contains basically the same information of the 
time dependence of the magnetization during relaxation, 
in particular, the higher is the sweep rate, the shorter 
is the effective "observing time" in the m(t) curve i 86 ! 87 
The results are shown in Fig. I10L along with the curve 
acquired by the SQUID and the values (full points) ex- 
tracted from the m(t) curves, obtained at different fields, 
considering t = 5000 sec. Regarding the measurement 
performed with the SQUID, one can assume SR = 
since any variation of H during the acquisition can be 
neglected. Furthermore, since the acquisition time is of 
some seconds, one has also to consider that the measured 
m value has already relaxed. A remarkable variation 
of the m(H) curves associated to different characteris- 
tic times is present, as a consequence of the relaxation 
in time of the magnetic moment. Furthermore, it is very 
interesting to note that the position of the second-peak 
moves towards low H values during the relaxation, which 
indicates that the magnetic field value at which the sec- 
ond peak in the J C (H) curve is obtained relaxes in time. 
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FIG. 10. Time evolution of the second peak position investi- 
gated by magnetization measurements performed by the VSM 
at different sweep rates (SR), by the SQUID as well as by 
magnetic relaxation measurements. 



In order to shed light on the mechanism that deter- 
mines the temperature and time dependence of the peak 
position, Hp ea k, we collected the couples of values (H pea k, 
Tripeak) which identify its location in the m(H) curves ob- 
tained at different temperatures by the SQUID and by 
the VSM and for different sweep rates, the results being 
shown in Fig. QTJ The inset at the bottom right of the 
same figure presents the calculated induction-field values 
(Bp ea k) corresponding to m pea k. It is very interesting to 
note that the pairs of values (H pea k, m-peak) - or equiva- 
lently (B pea k, m pea k) - obtained at different temperatures 
or sweep rates lie on the same curve and that, in some 
cases, the same pair of values (H pea k, m pea k) has been 
obtained at two different temperatures but for a differ- 
ent "observing time" during relaxation. The continuous 
line shown in the figure is the best fit curve obtained 
by assuming m pea k = a ■ ^oH^ eak . The best fit parame- 
ters, along with the corresponding statistical errors, are: 
a = 0.35 ± 0.01 emu/T; b = 1.31 ± 0.04. 

Since the induction field B is proportional to the num- 
ber of vortices present in the SC and the magnetic mo- 
ment identifies the critical current (and as a consequence 
the flux profile inside the sample), the pair of values 
(B pea k, mpeak) identify a specific configuration of the 
vortex structure. The fact that all the experimental data 
shown in Fig.[TT]fall in a unique curve suggests that there 
is a unique vortex configuration which determines the oc- 
currence of the peak effect at a given H pea k value, what- 
ever the temperature is in the range 4.2 K 8 K, cor- 
responding to 0.4 T c ~ 0.75 T c . In particular, the same 

eak ilTlpeak 

) pair can be found at a temperature T\ 
and observing time t\ or equivalently at a temperature 
T2 < T\ at a time £2 > t\- If one defines the vortex 
energy landscape as the mapping of all the possible con- 
figurations of the vortices in the sample, determined by 
B and m, and the corresponding energy level (E), the 
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FIG. 11. Pairs of values (H pea k,'m V eak) that identify the sec- 
ond peak position in the m(H) curve at different temperatures 
and characteristic times. Labels 1, 2, 3, indicate the measure- 
ments performed by the VSM with SR= 2T/min, lT/min 
and 0.05 T/min, respectively; label 4 indicates the measure- 
ments performed by the SQUID. The continuous line is the 
best fit curve obtained as described in the text. The inset at 
the bottom right shows the couples of values (B pea k,m vea k) 
corresponding to the peak effect. The inset at the top left 
is an imaginary drawing of the vortex energy landscape, the 
line indicating the relaxation path of the vortex configuration 
associated to the peak effect. 



observed behavior could be justified by supposing that 
there is a unique path in the vortex energy landscape 
that describes the relaxation of the vortex configuration 
associated to the peak effect and that the temperature 
does not particularly affect the energy landscape in the 
T range explored, at least in proximity of the path. The 
only effect of increasing the temperature is to allow the 
relaxation to start from a point in the path closer to the 
final equilibrium state. This result also indicates that the 
mechanism behind the second peak is related to a ther- 
mally driven process in the examined range of tempera- 
tures and fields. We are extending the present analysis of 
the peak effect properties to other superconductors, such 
as YBa2Cu307_,5, and the results will be discussed in a 
forthcoming paper. 



IV. CONCLUSIONS 

We have investigated the magnetic properties of 
FeTeo.7Seo.3 crystals in the superconducting state by 
magnetization and magnetic relaxation measurements. 
We have shown that pinning in FeTeo.7Seo.3 originates 
form spatial variation of the mean free path and is most 
likely related to the nanoscale chemical phase separation 
observed in Fei+^Tei-ySey on a scale of ~ 10 nm i 78 i 79 
Very interestingly, our results confirm that even if chem- 
ically inhomogenous Fei+^Tei-ySej, is electronically ho- 
mogeneous since pinning is not ruled by spatial fluctua- 
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tions of the critical temperature. From magnetic relax- 
ation measurements we have obtained indications that 
vortex motion develops in the single vortex limit even 
in magnetic field as high as 8 T, in agreement to what 
observed in the 122 Fe-based compounds. 4 iOn applying 
the magnetic field along the c axis, a clear peak effect 
in the m{H) curves has been observed, up to temper- 
atures near T c . The second-peak position varies with 
the temperature, following the same dependence as ob- 
served in YBa 2 Cu30y_5 Jfi The relaxation of the vortex 
configuration that determines the peak effect has also 
been studied. We have found that the pairs of values 
(H pea k, iTip ea k) that identify the second-peak position 
in the m(H) curves obtained at different temperatures 
and different relaxation times reconstruct a unique curve. 



This suggests that the vortex configuration that deter- 
mines the peak effect at a particular H value does not 
depend on the temperature in the range 0.4 T c 0.75 T c . 
It follows that the temperature influence on the vortex- 
vortex and vortex-defect interactions in proximity of the 
peak effect is negligible. 
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